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We present a first principles study of the series of multiferroic barium fluorides with the compo- 
sition BaA/F4, where M is Mn, Fe, Co, or Ni. We discuss trends in the structural, electronic, and 
magnetic properties, and we show that the ferroelectricity in these systems results from the "freez- 
ing in" of a single unstable polar phonon mode. In contrast to the case of the standard perovskite 
ferroelectrics, this structural distortion is not accompanied by charge transfer between cations and 
anions. Thus, the ferroelectric instability in the multiferroic barium fluorides arises solely due to 
size effects and the special geometrical constraints of the underlying crystal structure. 



I. INTRODUCTION 

Magnetoelectric multiferroics are currently attracting 
much attentioni^ These materials, which exhibit mag- 
netic and dielectric order in the same phase, can give 
rise to interesting coupling effects between the two fer- 
roic order parameters, with great potential for techno- 
logical applications. Currently, the main research effort 
in multiferroics is directed towards Bi-containing per- 
ovskites such as BiFe03 or BiMnOaji^ and towards 
both hexagonal and orthorhombic rare-earth manganites 
such as HoMn03 and TbMn03&i In the present article 
we revisit another class of magnetic ferroelectrics: the 
barium fluorides BaMF4, where M is Mn, Fe, Co, or Ni. 
This scries of compounds attracted considerable interest 
from the late 1960s until the early 1980s (see Ref. and 
references therein) but so far has not received much at- 
tention during the recent revival of multiferroic materials. 
We hope that our work will inspire further research in this 
interesting class of multiferroic materials, with the goal 
of exploring a broader range of non-oxide-based materials 
as candidates for magnetoelectric device applications. 

The barium fluorides BaMF 4 , M = Mn, Fe, Co, Ni, 
Zn, and Mg, form an isostructural family of compounds 
with the polar space group Cmc2\ (Refs. l9lHol) iS£ The 
corresponding crystal structure is shown in Fig. QJi. In 
this base-centered orthorhombic structure the transition 
metal cations are octahedrally surrounded by fluorine an- 
ions. Four of the six corners of the fluorine octahedra are 
shared with adjacent octahedra to form puckered sheets 
perpendicular to the orthorhombic b axis. These sheets of 
octahedra are separated by similar sheets of Ba cations. 
The structure is polar along the c direction, and ferro- 
electric switching has been demonstrated for M = Co, 
Ni, Zn, and Mg, but not for M = Mn and FeM The di- 
electric constants e c show an increase with temperature 
that is characteristic for a ferroelectric phase transition, 
but all crystals melt before a transition into the para- 
electric phase occurs. The melting temperatures range 
between 720-965 °C; the Curie- Weiss temperatures can 
be extrapolated from the temperature dependence of the 
dielectric constants and range between 810-1320 °C. 11 In 
BaMnF 4 an additional structural phase transition occurs 
at ^255 Ki— The resulting low-temperature structure is 



incommensurate along the c axis^ To facilitate a sys- 
tematic comparison between the different BaMF4 sys- 
tems we do not consider this low temperature structure 
of BaMnF4 in the present study and instead use perfect 
Cmc2i symmetry for all systems. 

In addition to the polar distortion, the systems with M 
= Mn, Fe, Co, and Ni exhibit antiferromagnetic order- 
ing below Tn ~ 50 - 120 °C4^ These materials therefore 
exhibit multiferroic behavior, which is the motivation for 
the present investigation. 

Here we present a comprehensive computational study 
of the structural, electronic, magnetic, and ferroelectric 
properties of BaMF4, with M = Mn, Fe, Co, and Ni, us- 
ing first principles electronic structure calculations. The 
goal is to elucidate the origin of ferroelectricity in these 
materials, and to understand trends, such as why the sys- 
tems with M = Mn and Fe do not exhibit ferroelectric 
switching. The origin of ferroelectricity is of fundamental 
interest, since BaMF 4 does not contain any ions, which 
are usually considered to be "ferroelectrically active", 
such as empty d shell cations or lone pair-active cations 
(e.g. Bi 3+ , Pb 2+ )pi£J£ We show that the ferroelectricity 
in these systems is due to the softening of a single polar 
phonon mode, which involves both rotational motions of 
the fluorine octahedra and polar displacements of the Ba 
cations. The instability is caused solely by size effects 
and geometrical constraints; no charge transfer between 
anions and cations occurs as a result of the structural dis- 
tortion. The BaA/F4 multiferroics therefore represent an 
example of proper "geometric ferroelectricity" , a mecha- 
nism that has been proposed as a possible way to incor- 
porate both magnetism and ferroelectricity in the same 
system^iiS Furthermore, our work represents the first 
ab initio study of a non-oxide multiferroic system, and 
we provide evidence that the LSDA+£7 method results 
in a good description of the electronic structure of these 
materials. 

This paper is organized as follows. We first describe 
the methods we use in our calculations, together with 
some technical details. We then present the results of 
our calculations for the structural, electronic, and mag- 
netic properties of all systems. Finally, we focus on the 
ferroelectric properties and analyze the mechanism un- 
derlying the polar structural distortions in these systems. 
We end with a discussion and summary. 
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II. COMPUTATIONAL METHOD 

All calculations in this work are performed using the 
projector augmented- wave method^ implemented in the 
Vienna Ab-Initio Simulation Package (VASP)i 20 i 21 
For the treatment of exchange and correlation we use 
both the local spin-density approximation (LSDA) — and 
the LSDA+[7 method in the formulation of Dudarev et 
al.fH which is equivalent to the standard form of Anisi- 
mov et aZ.f24 with the intra-atomic exchange parameter 
J set to zero (see Ref. I25J) . For the Hubbard parameter 
U e ff = U — J of the transition metal d states we use a 
typical value of 4 eV, and we test the sensitivity of our 
results with respect to the precise value of U c s where 
necessary (see also section UlIBfl . 

To obtain structural parameters we relax all ions un- 
til the Hellman-Feynman forces are less than 0.01 eV/A, 
and we adjust the lattice vectors such that all compo- 
nents of the stress tensor are smaller than 1 kbar. During 
the relaxations, except where otherwise noted, we impose 
the experimentally observed antiferromagnetic order (see 
Sec. IIII C|) . which doubles the size of the unit cell and 
reduces the space group symmetry to monoclinic. How- 
ever, no monoclinic distortion has been found experimen- 
tally, and in our calculations the corresponding effect, if 
present, is too small to be resolved unambiguously. We 
therefore neglect a possible monoclinic distortion and re- 
strict the lattice vectors to the orthorhombic symmetry. 
For test purposes, we also perform relaxations for differ- 
ent magnetic structures, e.g. ferromagnetic ordering. 

To obtain local densities of states we define spheres 
around the ions with radius 1.0 A for the M cations 
and 0.9 A for the fluorine anions. We use a 4x4x3 
Monkhorst-Pack k-point mesh (divisions with respect 
to the base-centered orthorhombic lattice vectors of the 
nonmagnetic unit cell) and a Gaussian smearing of 0.1 eV 
for Brillouin zone integrations. The plane-wave energy 
cutoff is set to 550 eV for relaxations and to 450 eV for 
all other calculations. To calculate the spontaneous po- 
larization we use the Berry-phase approach^SiSL 2 ^ and in- 
tegrate over 8 homogeneously distributed k-point strings 
parallel to the reciprocal c direction, each string contain- 
ing 8 k-points. For the calculation of the force-constant 
matrix we displace all the ions by 0.005 A, corresponding 
to symmetry-adapted modes. To exclude eventual non- 
linearities in the forces we repeat the calculation with 
displacements of 0.01 A, which leads to identical values 
for the force constant matrix. 



III. RESULTS AND DISCUSSION 

A. Structural properties 

Tablc|l]shows the structural parameters for all systems 
obtained here by using both LSD A and LSDA+[7 with 
U e g = 4 eV, together with experimental data. The overall 
agreement between calculated and experimental values is 



very good. The use of the LSDA leads to a typical un- 
derestimation of the lattice parameters between 1-4% 
whereas the internal structural parameters are very close 
to the experimental values. The LSDA+t/ method leads 
to a slightly larger equilibrium volume compared to the 
pure LSDA and therefore improves the agreement with 
experiment. This is very similar to what has been re- 
ported previously for several magnetic oxidesA^iM 

The calculated structural parameters shown in Table [I] 
are calculated for the experimentally observed antiferro- 
magnetic order (see Sec. IIII C|) . If a different magnetic 
structure is imposed during the relaxation (results not 
shown here), we notice small, but distinct, structural 
changes. For example, relaxation of BaNiF4 in a fer- 
romagnetic configuration leads to a 0.8 % increase in the 
lattice parameter b and up to 7 % change in the Wyckoff 
positions of those fluorine ions that mediate the mag- 
netic superexchange interactions in the M-F-M bonds. 
This indicates a certain degree of spin-lattice coupling in 
these systems which can give rise to phenomena such as 
spin-phonon coupling^ and magneto-capacitance, 



a) b) 




FIG. 1: (Color online) a) Projection of the BaAfF4 structure 
along the a axis. The M cations are octahedrally surrounded 
by fluorine anions, which form puckered sheets perpendicular 
to the b axis, separated by similar sheets of Ba cations. Adja- 
cent sheets are shifted relative to each other by half a lattice 
constant along the a direction, b) Corresponding centrosym- 
metric prototype structure. 

We also performed structural relaxations for all sys- 
tems within the corresponding centrosymmetric proto- 
type structure with space group symmetry Cmcm (see 
Ref. and Fig. QJ)), which can be obtained from the 
ground state structure by imposing an additional mirror 
symmetry perpendicular to the c axis. It has been sug- 
gested that the M-F(l)-M distance in this centrosym- 
metric Cmcm structure is too small to accommodate 
two M-F bonds, causing the structural instability in the 
BaAfF4 systemsj22 Here, F(l) is the fluorine anion con- 
necting neighboring octahedra along the c direction. The 
M-F(l)-M segment thus "bends" outward in order to in- 
crease the M-F(l) bond length, leading to the observed 
collective distortion of the octahedral network. This dis- 
tortion can be understood as collective alternate rota- 
tions of the fluorine octahedra around their respective 
centers with the rotation axis parallel to the a direction, 
accompanied by displacements of the Ba cations paral- 
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TABLE I: Structural parameters for Cmdl\ BaMF 4 , M—Mn, Fe, Co, Ni, obtained using the LSDA and the LSDA+C7 method 
with U e g=4 eV together with experimental data, a, b, c are the usual orthorhombic lattice constants, V is the corresponding 
volume, d represents the atomic displacements compared to the Cmcm centrosymmetric reference structure. All atomic 
positions correspond to Wyckoff-positions 4a: (0,y,z). 
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lei to the polar c axis. Table [I] shows the magnitude of 
the displacements d of both the Ba and F(l) ions leading 
from the centrosymmetric reference structure to the po- 
lar ground state structure. From these displacements it 
can be seen that the structural distortion increases over 
the series from M = Ni to M = Mn, which is also con- 
sistent with the experimental observations, and has been 
explained by a decrease in the ratio between the M-M 
distance and the minimum M-F bond length when M is 
changed from Ni to MniS 9 . 



B. Electronic structure 

The calculated total densities of states as well as the 
partial densities of the fluorine p and transition metal d 
states for all systems are shown in Figure |5] For M = 
Mn and Ni, even the use of the LSDA leads to an insu- 
lating solution. In the case of the Ni system the gap is 
small and is due to the crystal field splitting of the local- 
ized d states, whereas in the Mn system the gap is larger 
and is produced mainly by the strong exchange splitting 
between the two spin channels. The Fe and Co systems 
are metallic in LSDA, as expected for the d 6 and d 7 elec- 
tron configurations of the Fe 2+ and Co 2+ ions within the 
predominantly cubic crystal field. The very small band- 
width of the transition metal d states indicates an in- 
stability towards the formation of a Mott-Hubbard gap, 
and indeed the use of the LSDA+J7 method with U e s = 
4 eV leads to the formation of a large gap of about 2-3 eV 
for M = Fe and Co, and to a substantial increase in the 
width of the gap for the Mn and Ni systems. In all cases 
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FIG. 2: Total densities of states (gray shaded), partial fluo- 
rine p (shaded with thin diagonal lines) and transition metal 
d (black shaded) states (in states/eV), calculated within the 
LSDA (left panel) and by using U c s = 4 eV (right panel) for all 
BaMF4 systems (M = Mn, Fe, Co, Ni from top to bottom). 
Minority spin states are shown with negative sign. Zero en- 
ergy corresponds to the Fermi level (metallic systems) or the 
highest occupied state (insulating systems). 



the gap is between occupied and unoccupied d states of 
the transition metal ion. 

These results indicate that the use of the LSDA is in- 
adequate for the BaMF4 systems, whereas the LSDA+C7 
method with an appropriate U value leads to a good de- 
scription of the electronic structure. In this work, ex- 
cept where otherwise noted, we use t/eff=4eV, which 
is a typical value for transition metal cations in oxide 
materialsi^LiSSi Since the overall features of the tran- 
sition metal d states in the present fluorides are very 
similar to the oxide case, we expect that the same value 
of U is also appropriate for these systems. We point out 
that a variation of U within reasonable limits alters the 
electronic structure only qualitatively. Nevertheless, for 
quantities that are expected to depend critically on U , 
we always consider the explicit U dependence. 

It can be seen from Fig. [5] that all transition metal 
cations are in a high-spin configuration, where the local 
majority spin states are fully occupied, and the minority 
spin states are filled with 0, 1, 2, 3 electrons for M 
Mr. Fe, Co, Ni, respectively. In the case of the LSDA, 
the transition metal d states are energetically well sepa- 
rated from the fluorine p states, leading to only negligible 
hybridization between the two sets of states. The use of 
LSDA+J7 lowers the energy of the filled transition metal 
d states, leading to energetic overlap of these states with 
the fluorine p levels and a certain degree of hybridiza- 
tion, which is most notable for M = Ni. This increase 
in hybridization with increasing U could be an artifact 
of the LSDA+C7 method, which only corrects the tran- 
sition metal d states while leaving the fluorine p states 
unchanged. 

The densities of states for the centrosymmetric Cmcm 
structures (not shown here) are indistinguishable from 
those calculated for the ground state structures (both in 
LSDA and LSDA+L/). This indicates that the structural 
distortions in these systems do not lead to a significant 
change in covalent bonding, in contrast to the case of the 
ferroelectric perovskites, where the structural distortions 
lead to strong rehybridization between filled anion p and 
empty cation states42*^i We will further analyze the fer- 
roelectric instability of the BaMF4 systems in Sec. lIII Dl 



C. Magnetic properties 

The experimentally observed magnetic structure is 
shown in Fig. [3J 42 The magnetic M cations form quasi- 
two-dimensional puckered rectangular grids "parallel" to 
the a-c planes. Within each rectangular grid the mag- 
netic moments of nearest neighbors are aligned antipar- 
allel to each other. The coupling between adjacent grids 
effectively cancels, so that the magnetic order along the 
b direction is determined by the weak coupling between 
the next nearest neighbor planes that are ~14 A apart. 
This can lead to two different magnetic phases, in which 
the coupling along the b direction is either parallel (phase 
B) or antiparallel (phase A). Both phases have been ob- 




FIG. 3: (Color online) Magnetic structure of BaMF 4 , M = 
Mn, Fe, Ni (only magnetic ions are shown). For M = Co 
(phase A) the magnetic moments are oriented parallel to the 
c axis, but the relative orientations of the moments are the 
same. Gray lines outline the conventional orthorhombic unit 
cell; the puckered rectangular grids are shown in green. 




FIG. 4: Heisenberg nearest neighbor exchange coupling con- 
stants J a (circles) and J c (triangles). 



served for BaCoF4(^ and since the energy difference be- 
tween these two phases is very small, it is generally as- 
sumed that extrinsic effects such as defects etc. can lead 
to the preferred appearance of one or the other phase. In 
our calculations wc only consider phase A where the cou- 
pling between second nearest neighbor planes along the 
b direction is antiparallel. In all systems except BaCoF4 
the magnetic moments are oriented parallel to the b axis, 
whereas in BaCoF4 the magnetic moments are oriented 
parallel to the c axis (in both phase A and phase B)4£ 

The magnetic structure of BaMnF4 exhibits two dis- 
tinct variations of the structure described in the previ- 
ous paragraph and shown in Fig. [3J first, the antiferro- 
magnetic axis is rotated by 9 from the b towards the 
a direction^ and in addition all magnetic moments are 
slightly canted towards the c direction by about 0.1°, 
resulting in a very small magnetization^ 

We do not attempt to reconstruct the full magnetic 
structure of all BaAfF4 systems from first principles, 
but to check if the LSDA+C/ treatment of the electron- 
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FIG. 5: Heisenberg coupling constants for BaNiF4 as a func- 
tion of U e s (J a '- circles, J c : triangles). Open symbols cor- 
respond to calculations done with structural parameters ob- 
tained within LSDA. Filled symbols correspond to fully re- 
laxed calculations. 

electron interaction leads to a correct description of the 
magnetic properties of the multiferroic barium fluorides, 
we determine the nearest neighbor magnetic coupling 
constants for all systems. To do this, we first calculate 
the energy differences corresponding to different mag- 
netic configurations of the four spins in the unit cell. We 
then map the calculated energies on a simple Heisen- 
berg model, where we write the magnetic interaction as 
Eij — —2JijSi ■ Sj (si is the spin of cation i; J^ is the 
coupling constant between ions i and j) and consider 
only nearest neighbor interactions. Within our sign con- 
vention a negative Jij corresponds to antiferromagnetic 
coupling. 

Fig. 2] shows the calculated Heisenberg coupling con- 
stants J a and J c corresponding to pairs of nearest neigh- 
bor spins along the a and c direction, respectively, for M 
= Mn, Fe, Co, and Ni. All nearest neighbor couplings 
are antiferromagnetic, in agreement with the experimen- 
tally observed magnetic structure. The coupling becomes 
stronger from M = Mn to Ni, which is a consequence of 
the successive filling of the t<x g states (see Ref. l46|) . In ad- 
dition, the stronger hybridization between the transition 
metal d and fluorine p states in the case of the Ni sys- 
tem further increases the strength of the superexchange 
interaction, leading to particularly strong antiferromag- 
netic nearest neighbor coupling in BaNiF 4 . 

For BaMnF4 we obtain the exchange coupling con- 
stants J a — — 0.270 meV and J c = — 0.173 meV (using 
[/off = 4eV). This compares extremely well with exper- 
imental values that have been extracted from the mea- 
sured spin-wave dispersion: J a = — 0.282 meV and J c = 
— 0.197 meV (Ref. Il3j) . This very good agreement might 
be to some extent fortuitous, but can also be regarded as 
strong indication that the value of U e e indeed provides 
a good description of the electronic structure, since the 
magnetic superexchange coupling depends very strongly 
on U (see Fig. [5] and the following paragraph) . 

From the theory of superexchange it follows that 
the corresponding coupling constant is proportional to 
1/(7 One therefore expects a rather strong influence of 



FIG. 6: Energy differences AE per formula unit between the 
ferroelectric and the centrosymmetric prototype structures, 
calculated within the LSDA (circles) and by using C/ e ff=4eV 
(triangle). 



the Hubbard parameter on the strength of the magnetic 
coupling, if this coupling is mediated mainly by the su- 
perexchange interaction. Indeed, strong U dependence of 
the magnetic coupling constants has been found for vari- 
ous oxidesi 47 i 48 ' 49 Fig.Oshows the variation of J a and J c 
in BaNiF4 with the Hubbard parameter U c g used in the 
LSDA+J7 treatment of the electron-electron interaction. 
To see the pure electronic effect on the coupling constant, 
the structural parameters are kept fixed to their values 
obtained by the LSDA relaxation. It is apparent that the 
coupling strength significantly decreases with increasing 
[/ e fi, although the variation does not show an exact 1/U 
dependence. This is due to the fact that a variation of U 
in the bandstructure calculation also changes the over- 
lap of the wave-functions and therefore the transfer inte- 
grals, which are considered as constant in the theory of 
superexchange. 

The values corresponding to the relaxed structural pa- 
rameters for U e fi — 4eV are also shown in Fig. [3J It 
can be seen that there is only a negligible change of J a , 
whereas J c is reduced by 20 % due to the slightly different 
structure. This is due to the fact that the bond lengths 
and angles for the Ni-F(4)-Ni bond, which determines J a , 
are nearly identical in LSDA and LSDA+C/, whereas the 
Ni-F(l)-Ni bond angle, which determines J c , is reduced 
from 160.5° (LSDA) to 155.8° (U cS = 4eV). The Ni-F(l) 
distance increases slightly from 1.94 A (LSDA) to 1.96 A 
(C/ e ff = 4eV). 

In addition, if spin-orbit coupling is taken into account, 
there are small deviations from the collinear antiferro- 
magnetic structure^ which are due to the antisymmet- 
ric exchange or Dzyaloshinskii-Moriya interactional This 
small canting does not influence the structural and fer- 
roelectric properties studied in the present paper, and 
we therefore neglect all effects due to spin-orbit coupling 
here. 



D. Ferroelectric properties 

Fig.|S|shows the energy differences per formula unit be- 
tween the polar ground-state Cmc1\ structures and the 
corresponding centrosymmetric Cmcm structures, calcu- 
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TABLE II: Spontaneous polarizations (in /^C/cni 2 ) calculated 
with an ionic model using formal charges and by using the 
Berry-phase approach (U e g — 4 eV), together with available 
experimental data. 



ionic model 



Mn 
Fe 
Co 
Ni 



13.68 
10.34 
8.57 
6.34 



Berry phase 



13.60 
10.88 
9.02 
6.79 



Exp. [10] 



8.0 
6.7 



TABLE III: TO phonon frequencies u> of A a and B\ u symme- 
try for BaNiF4. x gives the contribution of the corresponding 
mode to the ground state distortion. 
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latcd within the LSDA and by using LSDA+C/ with U e g 
= 4 eV. The latter results in slightly larger energy differ- 
ences, but the overall trend is similar in the two cases. 
The energy difference is largest for M = Mn and smallest 
for M = Ni, consistent with the trend in the magnitude of 
the structural distortions (see Sec. IIII Af . The fact that 
no ferroelectric switching could be observed for BaAf F4 
with M — Mn and Fe 10 has been explained by suggesting 
a rather high energy for the intermediate centrosymmet- 
ric state in these two systems, which results from the dif- 
ferent equilibrium bond lengths compared to the systems 
with M = Co, Ni, Zn, and MgiSS This energy barrier is 
exactly the energy difference shown in Fig. EJ The values 
for M = Ni and M — Mn are comparable with the cor- 
responding energy differences between ferroelectric and 
centrosymmetric prototype structures in the perovskite 
ferroelectrics BaTi0 3 (« 18meV/f.u.) and PbTi0 3 (« 
200meV/f.u.), respectively^ Although it is not likely 
that ferroelectric switching actually takes place through 
the centrosymmetric reference structure, this energy dif- 
ference can be viewed as an upper limit for the switching 
barrier. It is apparent that a classification of the Co and 
Ni systems as ferroelectric and the Mn and Fe systems 
as merely pyroelectric is probably not very instructive, 
since there is only a gradual quantitative difference be- 
tween these systems. It seems feasible that ferroelectric 
switching could be achieved at least for the Fe system if 
highly resistive samples were available that could sustain 
higher electric fields. 

From the displacements of the ions between the cen- 
trosymmetric prototype and the polar ground state struc- 
ture it is possible to calculate an estimate of the spon- 
taneous polarization using the formal charges of the ions 
(2+ for the cations and 1- for the fluorine anions). Ta- 
ble[n]lists the polarization values calculated from such an 
ionic model together with the corresponding values ob- 
tained by ab initio calculation using the Berry-phase the- 
ory of electric polarizationi2§i2L2& The good agreement 
between the Berry-phase result and the ionic model in- 
dicates that the Born effective charges^ are very close 
to their formal values for all systems. This is consis- 
tent with the lack of significant changes in the electronic 
densities of states between the centrosymmetric and dis- 
torted structures as described in Sec. IIII Bl Thus, no 
significant charge transfer occurs as a result of the polar 



distortion. This again indicates that the ferroelectricity 
in these systems is of different origin than in the pro- 
totypical perovskite ferroelectrics where charge transfer 
between the transition metal d and the anion p states is 
crucial for stabilizing the ferroelectric stated 

To obtain further insight into the ferroelectric insta- 
bility of these systems, we calculate the TO phonon fre- 
quencies of the centrosymmetric structure for BaNiF4. 
Since we are only interested in the phonon modes that 
are compatible with the observed ground state symme- 
try (only these modes can give rise to the correspond- 
ing structural distortions), we transform the dynamical 
matrix to block-diagonal form using symmetry-adapted 
modes, and diagonalize only the blocks corresponding to 
the irreducible representations A g and B\ u of space group 
Cmcm. The nonpolar A g phonons do not change the 
symmetry of the system, whereas the infrared active B\ u 
modes reduce the symmetry to the ground state Cmc2i 
space group. There are 5 A g and 7 B\ u modes. The cal- 
culated phonon frequencies are listed in Table ITTT1 One of 
the 7 B\ u modes corresponds to an acoustic mode with 
zero frequency at the T point and is not included in Ta- 
ble EH 

We find that there is one unstable B\ u mode with 
lu — i-BScm^ 1 . If we decompose the structural distor- 
tion leading from the centrosymmetric reference struc- 
ture to the relaxed ground state structure into contribu- 
tions from the different phonon modes, we see that the 
unstable B\ u mode is responsible for about 87 % of the 
final distortion. This indicates that the ferroelectricity 
in the barium fluorides originates from the softening of 
a single polar phonon mode which "freezes in" to form 
the ferroelectric ground state of the system. In contrast 
to the well-known case of perovskite ferroelectrics such 
as BaTiOs or PbTiOa no charge transfer between an- 
ions and cations occurs as a result of the displacements, 
which would lead to anomalous values of the Born effec- 
tive charges, but instead the structural distortion in the 
BaAfF4 systems is driven purely by size effects, which, to- 
gether with the special geometric connectivity realized in 
these compounds, leads to an inversion symmetry break- 
ing and the appearance of a spontaneous electric polar- 
ization. 
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IV. DISCUSSION AND SUMMARY 

It is well known that size-effects can drive structural 
distortions. For example in the perovskite system the 
relative ionic radii of the A and B cations determine the 
stability towards tilts of the anion octahedra. 53 

The idea that ionic polarizability is, in general, not 
necessary for achieving ferroelectricity was already for- 
mulated in Ref. |?3 when a ferroelectric ground state 
was predicted for perovskite-like NaCaF3 on the ba- 
sis of atomistic simulations using Gordon-Kim pair- 
potentials^ Since the possibility of covalent bond for- 
mation is excluded within this method, it became clear 
that a different mechanism is driving the ferroelectric dis- 
tortion in NaCaF 3 . Although the underlying mechanism 
was not investigated any further, it was assumed that 
ionic polarizability would significantly enhance the ferro- 
electricity in NaCaF3. The present study shows that this 
does not necessarily have to be the case. Ferroelectricity 
was also predicted recently for a variety of perovskite-like 
fluorides, such as LiMgF 3 , LiNiF 3 , and NaCdF 3 ^ilt is 
reasonable to assume that the ferroelectricity in all these 
fluoride systems is driven purely by size effects and that 
no charge transfer between anions and cations occurs. 

It was pointed out recently that such a mechanism for 
ferroelectricity is compatible with the simultaneous oc- 
currence of magnetic ordering^ Using first principles 
electronic structure calculations, the hexagonal multi- 
ferroic YMn0 3 was identified as an example for such 
a scenario. In this material no charge transfer occurs 
as a result of the structural distortion, and the dynami- 
cal charges are very close to their formal values ^1 This 
suggests that in YMn0 3 , similar to the BaMF4 systems 
discussed in the present paper, ionic polarizability does 
not have any effect on the ferroelectric properties. A 



subsequent analysis of the phonon modes of the corre- 
sponding centrosymmetric reference structure suggested 
that YMn0 3 is in fact an improper ferroelectric, where 
the ferroelectricity is due to the symmetry-allowed cou- 
pling of a stable polar zone-center phonon mode to an 
unstable zone-boundary modeS 

The present study shows that such "geometric ferro- 
electricity" is not restricted to improper ferroelectrics, 
but that size effects can also directly lead to instable po- 
lar zone center phonon modes. Thus, BaNiF4, and with 
it the whole class of BaMF4 multiferroics, represent the 
first confirmed example of a proper geometric ferroelec- 
tric. 

In summary, we have presented an ab initio study of 
the structural, electronic, and magnetic properties of the 
multiferroic barium fluorides BaMF4 with M = Mn, Fc, 
Co, and Ni. The ferroelectricity in these systems is due 
to a single unstable polar zone-center phonon. The in- 
stability is triggered solely by size effects, and no charge 
transfer occurs as a result of the structural distortion. 
In addition, we have shown that the LSDA+L 7 method 
results in a good description of the electronic structure 
of these systems. Our work represents the first ab initio 
demonstration of a proper geometric ferroelectric and the 
first such study of a non-oxide multiferroic. 
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